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GUT-scale parameters assuming universal boundary conditions. We give pre- 
dictions in this model for the inclusive 6 — > 57 branching fraction and in- 
vestigate the impact of non- universal scalar mass boundary conditions. We 
find our results do not depend significantly on the value of the GUT-scale 
trilinear coupling A . The small tan/3 region favors predictions for the in- 
clusive 6 ^ 57 branching fraction close to that of the Standard Model value. 
Nevertheless forthcoming experimental results can eliminate some regions of 
the GUT parameter space. 
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It is now widely appreciated that the b ^ s'y inclusive branching fraction is potentially 
sensitive to new physics beyond the Standard Model (SM), such as charged Higgs bosons 
and supersymmetric particles [|1]-^. The inclusion of QCD corrections to the one- loop ampli- 
tudes significantly enhances the decay rate and these calculations have been steadily refined 
l^-j^. The SM amplitude is dominated by the loop with an intermediate weak boson W 
and top quark t, and the branching fraction depends on the mass of t for which there is now 
some direct evidence |^. With supersymmetry (SUSY) there are additional contributions, 
the most important coming from loops with charged Higgs if ^ and t, charginos Xi cind 
stop tj. The phenomenology can be complex in general, because of the potential for can- 
cellations between different loop contributions, and because of the large number of possible 
parameters. However, considerations of supersymmetric grand unification (GUT) provide 
great simplification, with many low-energy parameters determined by just a few GUT-scale 
parameters through the renormalization group equations. In this Letter we investigate the 
predictions for b —>■ s'-f based on the ROE analysis of Ref. [0], and show that measurements 
of the B -^ Xs7 inclusive decay rate can exclude otherwise allowed regions for the GUT 
parameter space. 



An attractive possibility is the realization of an infrared fixed-point solution W^ for the 



top-quark Yukawa coupling Af in supersymmetric GUT models |]T0|,|T2|-[T^ . In this scenario 



the top-quark mass rrit is related to the ratio of the two Higgs vacuum expectation values 
'y2/'f^i = tan/3 by ||12 



mt(pole) - (200 GeV) sin/3 . (1) 

In the following discussion we use the value mt{pole) = 168 GeV and as(M^) = 0.120, for 
which the RGE solutions for the sparticle mass spectra were constructed in Ref. |T^; this 
is consistent with the value rrit = 174 ± IOIJ2 experimentally observed from the top-quark 
candidate events p|. The universal soft SUSY-breaking parameters at the GUT scale Mq 
are uiq, mi/2, /i, A, B; here itiq is the scalar mass, mi/2 is the gaugino mass, /i is the 
Higgs mass mixing, and A, B are trilinear, bilinear scalar couplings respectively. In the 



ambidextrous approach of Ref. |T0| (see also Ref. [|T5|) used to solve the RGE integration, 
mi/2, ^^0, and A are input at the GUT scale. Near the fixed point, the value of the trilinear 
coupling At that enters into squark mixing is drawn to a fixed value (for a given a^lM^)) 
|1^,|T6|, so that the dependence of the rate on the initial value A'-' is minimal. These soft 
SUSY-breaking parameters are evolved down to the electroweak scale M^, and the values 
\fi{M^)\ and B{Mz) are determined from the requirement that the electroweak symmetry 
be broken radiatively. By this procedure, all sparticle masses and couplings are determined 
in terms of GUT parameters mo and mi/2 for a specified sign of /i. Thus we can predict 
the b —>■ s'-f decay rate over the region of (mo, rni/2) parameter space allowed by other 
constraints, which are experimental lower bounds on sparticle masses and a naturalness 



criterion taken to be |/i| < 500 GeV in Ref. |1Q 



For calculating the 6 — i> 57 decay rate, we use the latest analysis of QCD corrections by 
Buras et al. 0, which evaluates the full 8x8 anomalous dimension matrix. We do not take 
into account the small effect of running the scale of the decay process from Q = rrit down 
to Mw (estimated in Ref. |T^), beyond the usual effect of running from Q = Mw down to 
rri},. The QCD effects also receive theoretical corrections from other large mass splittings in 
the super symmetric spectrum |T^ . 

The ratio oi r{b —^ 57) to the inclusive semileptonic decay width is then given by 



r(6^s7) 6a \VZVt 
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where a is the electromagnetic coupling. The phase-space factor p and the QCD correction 



factor A for the semileptonic process are given by p = 1 — 8r^ + 8r^ 



24r*^ ln(r) with 



r = mjmb = 0.316 ± 0.013 and A = 1 - f /(r, 0, 0)a^{mb)/n with /(r, 0, 0) = 2.41 ^ 
The formulas for the renormalization group coefficient in Eq. (|^) are 0,^ 
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with Zi = {ni/ni^^y and Z2 = {niti,/ni^_.y and the matrices T, U and V diagonahze the top 
squark and chargino mass matrices (see Refs. [0,|l^). The parameter m is the mass of the 
(degenerate) squarks from the first two generations. 

The coefficients from the 8x8 anomalous dimension matrix are [161 
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hi 
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i' i' -i' 0.4086, -0.4230, -0.8994, 0.1456) 
^^2^^ 3 -^, -0.6494, -0.0380, -0.0186, -0.0057) 
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The formulas for the loop contributions f^^g are given in Ref. [^. 

The neutralino and gluino loop contributions have not been included in the calculation 
of the 6 — i> S7 inclusive rate. However, these contributions are known to be small, especially 
in the low tan/3 regime. 

Our results are displayed in Figs. |l|-0. Figure |l] and Figure ^ show the predicted value 
of the inclusive branching fraction B{b —* S7), that is equal in the spectator approximation 
to the physical meson branching fraction B{B ^ -^s7); the upper part of these figures 
gives the branching fraction versus m^ for fixed niii2 , while the lower part gives curves 
versus niii2 for fixed ni^ . The same information is displayed in Figure ^, in the form of 
contours of constant B(h -^ 57) in the {nii/2,nio) parameter plane. The region is bounded 
by experimental limits on the chargino mass m ± and the lighter stop mass m^^ and the 



lighter higgs mass rrih together with the naturalness boundary \fi\ = 500 GeV. The shaded 
region is allowed by these constraints. 

It has come to light in recent months that the leading order QCD calculations are some- 
what inadequate for comparison with the current experimental data; hence a brief discussion 
of the theoretical uncertainties in warranted. The (rrife)^ dependence of the rate for b ^ s'j 
is customarily removed by taking the ratio with semileptonic decays. This, however, does 
not entirely eliminate the theoretical uncertainty due to rrih because the 6-mass still enters 
into the phase space factor which depends on the ratio mb/mc- The ratio of CKM matrix 
entries is taken to be the central value [^ 






0.95 ± 0.04 . (9) 
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There is even more substantial uncertainty in the higher order (three-loop) contributions. 
It has been argued [^ that the unknown next-to-leading QCD corrections yield a ±25% 
uncertainty in the branching fraction. This number is arrived at by varying the unphysical 
scale /i from ?7if,/2 to 2?7ifc, and is particularly large since the process 6 -^ 57 is dominated 
by the QCD corrections. The calculation of the next-to-leading corrections is a formidable 
task, and progress in reducing this theoretical uncertainty is not anticipated in the near 
future. Furthermore, the value of ^^(M^) is not known very precisely, and for ^^(M^) in 
the range 0.12 ± 0.01, the b ^ sj rate changes by roughly 10%. Finally, the experimental 
error in the measurement of the semileptonic decay rate must be considered. 

Most analyses of the 6 ^ 57 process have focused on the case where the scalar masses 
are universal at the GUT scale. It has recently become fashionable to consider situations 
where the scalar masses take a more general form but still satisfy the stringent bounds on 
the flavor changing neutral currents of the first two generations. This nonuniversality might 



be generic |2^ , or it might arise from the running of the scalar masses from a universal value 
near the Planck scale to the GUT scale [^ . When the gauge group is broken to one of lesser 
rank and non-universality in the soft-supersymmetry breaking terms exists, then there are 
in general additional D-term contributions to the scalar masses p3| , p4| . One must be careful 



to include the often-neglected contributions to the renormalization group equations of the 



soft-supersymmetry breaking parameters that arise when non-universality is assumed p5 
The combination 

5 = m^^ - ml^ + Tr[M^^ - Mi^ - 2M^^^ + M^^ + M|J , (10) 



satisfies the (one-loop) scaling equation |g4 



so that if it is zero at some scale, for example the GUT scale, then it is zero for all scales. 
The renormalization group coefficient hi = 33/5 in the minimal supersymmetric model. We 
have used the two-loop renormalization group equations for which three groups are now in 
complete agreement |2^. The supersymmetric spectrum that results will be different from 
that obtained with universal masses at the GUT scale. In a particular example based on 
the minimal SU{5) GUT investigated recently by Polonsky and Pomarol ||22|, each of the 
sfermion (10 and 5) and the Higgs masses is separated at the GUT scale due to the evolution 
from the Planck scale. If one takes the sfermion mass to be twice that of the Higgs scalar 
mass at the GUT scale, the approach of B{b -^ sj) to the Standard Model value is faster 
than in the universal case as one increases the scalar masses: see Figures 4 and 5. This 
results from the contributions of the comparatively heavier top squarks which enter into the 
chargino loop diagrams. 

One can also consider the nonuniversality of the two Higgs doublet mass parameters 



| 27| . In this case one finds that the value of /x needs to be adjusted to achieve the correct 



electroweak symmetry breaking as described by the (tree-level) expression 

1 r, m?TT — m?TT tan^ d r, , , 



The charged-Higgs mass can be increased or decreased, but since /i is large for small tan/3, 
the mass is guaranteed to be large as well. Consequently the impact on the b ^ s^y rate is 
small. In some small regions of parameter space, the mixing between the top squarks can 
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be enhanced. However the overall qualitative behaviour depicted in Figures 1 and 2 is not 
changed substantially. 

The effects of the nonuniversality of the symmetry breaking terms on B{b -^ sj) in 
the low tan/5 region can be summarized as follows: (1) There is little dependence on the 
chargino mass and couplings in the low-tan f3 fixed-point region since fi is large and dominates 
the chargino mass matrix. Therefore the mass of the lightest chargino is approximately 
M2 = ^mi/2, and its coupling is predominantly gaugino. (2) Consequently, the dominant 
effect comes from the top squark masses and mixing. Typically increasing ttiq relative to 
rriQ increases the top squark mass and pushes the 6 — > S7 rate closer to the Standard Model 
value. (3) The dependence on A'^ is minimal near the top Yukawa coupling fixed point 
because the low energy value At is driven to a fixed value |lT^ , p!6| . It has pronounced effects 



only in regions where there is large top squark mixing. This results in the lightest top 
squark mass being suppressed (below the experimental limit and one is even in danger of 
the mass-squared going negative) and an enhancement of the chargino contribution to the 
6 ^ 57 inclusive rate. 

To summarize, we have evaluated the inclusive decay branching fraction B{b ^ 57) 
through the allowed regions of SUSY-GUT parameters, for solutions with small tan/3 near 
the At fixed point, following Ref. [|lOl- Our results show that 

(a) In the low tan/3 fixed point region, the inclusive rate for b —>■ s'j is typically close to that 
expected in the Standard Model. It is larger (smaller) than the Standard Model value for 
yU > (/i < 0) when the difference is substantial. 

(b) An accurate measurement of B{b — > 57) must be coupled with reduced theoretical 
uncertainties in order to constrain the parameters mi/2 and mo, or to test the universality 
of the soft-supersymmetry breaking parameters. 

(c) The contours in Figure 3 illustrate the kind of constraints that may be expected. 
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FIGURES 

1. Predicted values oi B{h -^ sj) (a) versus ttiq for fixed values of mi/2 and (b) versus 
mi/2 for fixed values of ttiq , in the low tan f3 SUSY-GUT solutions of Ref . [jlO[ with 
^ > 0. The curves are shown in intervals of 10 GeV in mo and mi/2; each curve in (a) 
corresponds to a vertical line in (b). 

2. Predicted values of B{b -^ 57) (a) versus uiq for fixed values of mi/2 and (b) versus 
mi/2 for fixed values of ttiq , in the low tan (3 SUSY-GUT solutions of Ref. []TU[ with 



/i < 0. The cross indicates the point at which the lighter top squark mass-squared 
becomes negative. 

3. Contours of fixed B{b -^ 57) in the (^«i/2,'"^o) parameter plane, for (a) yU > and 
(b) fi < 0. The shaded areas are allowed by direct experimental and naturalness 



constraints 10 



4. Predicted values oi B{h ^ 57) for yU > (a) versus mg = 2m^ for fixed values of mi/2 
and (b) versus Tni/2 for fixed values of m^ . 

5. Predicted values of 5(6 ^ 37) for yU < (a) versus ml = 2m^ for fixed values of mi/2 
and (b) versus mi/2 for fixed values of m^ . 
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